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Abstract. We provide an analytical approximation of a periodic solution of the three 
body problem in celestial mechanics, the so-called figure eight solution, discovered 1993 
by C. Moore. This approximation has the form of a Fourier series whose components 
are in turn Taylor series w. r. t. some parameter. The method is first illustrated by 
application to two other problems, (1) the problem of oscillations of a particle in a cubic 
potential that has a well-known analytic solution in terms of elliptic functions and (2) 
periodic solutions and corresponding eigenvalues of a generalized Mathieu equation 
that cannot be solved analytically. When applied to the three body problem it turns 
out that the Fourier-Taylor series, evaluated up to 30th order, represents un-physical 
solutions except for a particular value of the series parameter. For this value the series 
approximates the numerical solution known from the literature up to a relative error 
of 1.6 X 10-3. 


1. Introduction 

Analytical approximations provide an access to problems that cannot be solved 
analytically and have some advantages in comparison to pnrely nnmerical solntions. 
Let us consider periodic solutions x(t) of an n-dimensional differential equation of the 
form 

X = f(x) . (1) 

An obvious attempt to solve ([1]) is by means of Fourier series for x(t). However, if 
f is non-linear one encounters the problem to evaluate the r. h. s. of ([T]). Already 
the quadratic term of a Taylor expansion of f(x) involves the convolution of two inhnite 
series of Fourier coefficients. Fortunately, in physical problems it often happens that the 
Fourier coefficients of the pertinent quantities decay exponentially with order. Hence one 
could try to truncate the inhnite Fourier series thereby avoiding the inhnite convolution 
problem. In order to get a better approximation one has to increase the length of the 
hnite Fourier series stepwise and one would like to have an iterative solution algorithm 
that uses the previous results for the coefficients at each step. The method of Fourier- 
Taylor (FT) series is a systematic approach that satishes these requirements. Moreover, 
it solves another problem, namely that the period of the oscillation that has to be used 
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in the Fourier series ansatz is often unknown and can only be numerically calculated. 
For the above problem the general FT series would read 


OO 


exp{inu:t) 

nGlj m=\n\ 

OO 

(2) 

m=0 

(3) 


Hence each Fourier component of x(f) of order n is a Taylor series w. r. t. a parameter 
A that starts with a term proportional to A”. Put differently, x(f) is written as a Taylor 
series w. r. t. A such that each Taylor coefficient of A"* is a hnite Fourier series at most of 
order m. Hence, if the r. h. s. of (II]) is expanded into a Taylor series w. r. t. A we only en¬ 
counter convolutions of hnite order. A comparison of the coefficients of A™ expfinut), 
n = 0,..., m yields a sequence of equations that can be used to iteratively determine 
the unknowns Anm and Um- This comparison is performed by expanding u only at the 
1. h. s. of dll); at the r. h. s. of dlj) the oj occurring in expfinut) is left intact. 

Up to now, A appeared as a purely formal parameter that is only used for book-keeping. 
However, this turns out to be problematic. Nothing prevents one to use, say, sinh(A) 
instead of A as the series parameter. This transformation would also inhuence the un¬ 
known coefficients and hence these coefficients a priori cannot be uniquely determined. 
Already counting the number of equations and unknowns indicates this kind of under¬ 
determination. The way out of this problem is a “concretization” of A, i. e. , A has to be 
chosen as a parameter with a concrete meaning, e. g. the amplitude of the oscillation. 
This will modify the FT ansatz (I2|) and, hopefully, give unique solutions for A^n and 
Ujn- Nevertheless, the hrst coefficients are often ambiguous and some choice has to be 
made, see the examples of the following sections. This is due to the effect that the very 
hrst equations are quadratic or of higher order and become linear ones only after a few 
steps. 

The method of FT series mostly cannot be successfully applied without carefully 
considering the proper choice of the variables and of the differential equation. Often one 
has some information about the solutions that can be used to simplify the ansatz (I2|) and 
(1^ . This is also exemplihed in the following sections. Generally speaking, the method 
of FT series can be viewed as a generalization of the method of linearizing equations 
like dl]) and, similarly, its quality of approximation depends on the size of the parameter 
A and of the maximal order of truncation. Since the basic idea is very simple we would 
not be surprised to learn that the method has already been used before, although we do 
not know of any application in the literature. One reason for this may be that the use of 
high truncation orders requires computer algebra software that has been only available 
during the last decades. To our best knowledge, the method has hrst been used in a 
lecture of one of the authors on non-linear wave equations [1]. 
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The paper is organized as follows. In order to illustrate the concrete application of 
the FT series method and to demonstrate its performance for a simple but non-trivial 
example we consider, in section [2l the oscillations of a particle in a cubic potential (an- 
harmonic oscillator). Since the analytical solution of this problem is well-known it can 
be utilized to test the quality of approximation provided by the method. We expand the 
FT series up to 24th order and obtain a precision of 5 x 10“^® for a medium amplitude 
that is 1/6 of the maximal value obtained for the aperiodic limit case. 

Secondly, in section |3] we consider periodic solutions of a generalized Mathieu equation of 
Hill’s type. Possible physical applications include the Schrodinger equation of a particle 
in a periodic potential and diffraction of electromagnetic waves in a periodic grating. 
In contrast to the an-harmonic oscillator this is a linear problem such that the Fourier 
series ansatz leads to an inhnite matrix problem. Periodic solutions exist only for certain 
values of a parameter a usually called “characteristic values” of the differential equation. 
We will refer to these as “eigenvalues” in accordance with the nomenclature commonly 
used in physics. Interestingly, the FT series leads to some kind of perturbation series for 
the eigenvalues and eigenfunctions without using the methods of perturbation theory. 
For a similar approach to holographic diffraction see [2] . 

As the main application we consider, in section HJ the periodic solution of the three body 
problem that has been found two decades ago [3] and mathematically investigated in [1] , 
but where an analytical solution is not available. Three equal masses perform a periodic 
motion following the same orbit that has the form of the hgure eight with a constant 
relative time delay of T/3 (“choreography”). Without loss of generality the period T 
can be chosen as T = 27r. Usually the FT series method yields a one-parameter family 
of periodic solutions as in the example considered in section |2l For the hgure eight 
solution we also obtain such a family but it will be completely un-physical due to its 
violation of angular momentum conservation, except for a particular value Ai ~ 2.07552 
of the parameter A. For this value the hgure eight solution numerically calculated in |1] 
is reproduced with a relative error of less than 1.6 x 10“^ if we expand the FT series up 
to 30th order. We close with a summary and outlook. 

The detailed results of the FT series method for the three examples are too complicated 
to be presented in a paper and will be given in three MATHEMATICA 10.1 notebooks 
that can be downloaded from [5]. Further applications of the FT method will be given 
elsewhere [6]. 

2. FT series for oscillations in a cubic potential 

We consider the oscillations of a particle with unit mass in a potential 

V{x) = 2x {3a^ + A^) - 2x^ , (4) 

and total energy 

E = —4a(a — A)(a -|- A) , 


(5) 


FT series for the figure eight solution 


4 


Energy 



Figure 1. Illustration of the potential V{x) according to ([4]) and the values of a = 1 
and A = 1/2 corresponding to a total energy oi E = —3. The particle performs 
non-linear oscillations between the positions Xmin = —1.5 and Xmax = —0.5 . 


X 



Figure 2. Oscillation of a particle in the cubic potential of figure [T] according to the 
analytical solution (dashed curve) and the FT series of 24th order (red curve). 


where a and A are positive parameters. The potential is shown in hgure[I]for the choice 
of a = 1 and A = 1/2. 

Energy conservation yields the differential equation 

= 2{E — V(x)) = 4(x — 2a){x — A -b a)(x + A -b a) (6) 

= Ax^ - g 2 X - g 3 = P{x) . (7) 

Note that this differential equation differs from the form given in ([1]), but it has the 
advantage that its r. h. s. is a polynomial in x and hence has a very simple Taylor 
expansion. 

It is easily shown that P{x) is positive in the range Xmin = —a — X < x < —a + X = x^ax 
and hence the particle performs oscillations between the positions x^in and Xmax- In 
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Figure 3. Difference between the two solutions displayed in figured] 


the limit A —?• 0 we obtain the asymptotic solution 

x{t) = — a + X cos ojt + (8) 

a; = 2\/^ + 0{X‘^) . (9) 

Here and in what follows we chose the point of zero time such that the maximal ampli¬ 
tude is assumed at f = 0. 


For the general solution we make the following FT series ansatz: 

OO OO 

x(t) = AnmX^ COS U CO t 

n=0 m=n 

OO 

U =J2^mX^. 

m=0 


( 10 ) 

( 11 ) 


The FT series contains only cos-terms since x{t) is an even function for the above 
choice of the zero time point. This ansatz is inserted into the differential equation ([6]). 
Both sides of ([6]) must have the same coefficients of the terms A™ cos nut. This yields a 
sequence of equations that are used to determine the unknowns Um and An m to arbitrary 
order that is only limited by computer resources. It is advisable to do the first steps 
“by hand”. 

The coefficient of A° of the 1. h. s. of ([6]) vanishes, whereas the r. h. s. of ([6]) gives 
4(Hoo —2a)(Hoo + a)^. We choose the solution Aqq = —a in view of (IH]). The coefficient of 
A^ vanishes at both sides of ([6]). Comparison of coefficients of A^ cos nut for n = 0,1,2 
yields the following three equations: 

0 = 12a{2Al, + Al,-2) + Al^ul, ( 12 ) 

0 = a Aqi All , (13) 

0 = Aii{ul-I2a) . (14) 


We choose the solution Hqi = 0, An = 1, Uq = 2\/ofi in accordance with ([8]) and ([9]). 
After these choices the comparison of coefficients of orders m > 3 yields unique results 
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at least up to order 24. It turns out that the coefficients of the powers A™' are only 
non-zero if m increases in steps of 2 in accordance with the above solution Aqi = 0. 
The detailed results are too complicated to be presented here and can be found in [5]. 
We only give the next order corrections to ([8]) and ([9]): 


y2 

x{t) = — a + ——h ( A 


12a 


A3 A 

I cos ut 


192a2 


A2 

12a 


ui 


= 2\/^ 


A2 


+ 0{\^) . 


cos 2 uj t + 0{X^), (15) 

(16) 


8\/3a3/2 

Finally we compare our FT series solution of order 24 with the analytical solution 
for the special values a = 1 and A = 1/2. For the frequency u and the corresponding 
oscillation period T we obtain 

27r 

3.4458763864722597, T = — ^ 1.8233925429960192 . (17) 


OJ 


ui 


The period can also be obtained by the following integral 

1 

2 1 


T = 2 


= 2 


dx 


^4(a;-2) {x + ^) (x + |) 

= 1.823392542996019014655... , 




(18) 


(19) 


where K{m) denotes the complete elliptic integral of hrst kind, see Ch. 17. We note 
that the relative error is of order lO”^^. 

Upon defining 


Ml = 


4(a; - 2) (x |) [x + |) 
and evaluating the coefficients of P{x) 

5^2 = 4 (3a^ -|- A^) = 13, gs = S (a^ — aA^) = 6 , 
we obtain after some standard transformations 


7 V 7 


dx = -. -K{-\+z. -K{- \ ,(20) 


7 V 7 


( 21 ) 


x{t) =V{t - ui,g2,g3) , ( 22 ) 

where V{z] g 2 , gfi) denotes the Weierstrass elliptic function, see [7], Ch. 18. The 
coincidence of this analytical result with the above FT series of order 24 is very close, 
see the hgures [2] and [3l Of course, the quality of the approximation would decrease 
when approaching the aperiodic solution corresponding to A = 3a. We will not closer 
investigate this problem since our intention was only to provide an example where the 
method of FT series gives reasonable results. 
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Figure 4. The lowest eigenvalue a for odd solutions of the generalized 
Mathieu equation (1231) as a function of the parameters A and /i. For /r = 
0 the lowest eigenvalue can be calculated by the MATHEMATICA command 
MathieuCharacteristicB[l, A^/2], see the blue curve. 





Figure 5. A special periodic solution fix) of the generalized Mathieu equation 
(1^ with X = 2, ^ = —0.3 according to the FT series (blue dashed curve) and the 
corresponding numerical solution (red curve). 
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Figure 6. The difference 6 ip between the two solutions show in figure [51 

3. FT series for the generalized Mathieu equation 

As a further application of the FT series approximation we consider a linear differential 
equation that has periodic solutions with period (or wavelength) T = 2tt for the 
eigenvalues a, namely a generalized Mathieu equation of the form 

(F 

= —a + cos(2a;) + A"^/i cos(4a;) . (23) 

dx^ 

We have called the independent variable x and the unknown function 'ip{x) because of the 
association with the one-dimensional Schrodinger equation for a particle in a periodic 
potential. For /i = 0 this is the ordinary Mathieu equation [7], Ch. 20, with = 2q 
and its periodic solutions including the corresponding eigenvalues are known special 
functions that be calculated by using computer-algebraic commands. However, the 
corresponding eigenvalues and solutions of the generalized Mathieu equation could only 
be numerically calculated, e. g., by truncating an inhnite-dimensional matrix problem. 
Since we only want to illustrate the application of the FT method it will suffice to 
consider odd solutions corresponding to the lowest eigenvalue of (|23ll . For A = 0 this is 
the solution 'ip{x) = A sinx with the eigenvalue a = 1, where we have set the amplitude 
to A according to the FT ansatz (j2]). In the same way the pre-factors A^ and A^ in fl23l) 
are dictated by the FT ansatz, whereas /r is a free real parameter. The FT ansatz for 
odd solutions of fl2^ is hence chosen as 

^/’(x) = A sin(x)-b ^ ^ A™ sin(ux), (24) 

n=3,5,... m=n,n+2,... 

= 

m=0,2,4,... 


a 


(25) 
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Figure 7. Illustration of the periodic figure eight solution of the three body problem. 
Without loss of generality the period is chosen as T = 27r. At time t = 0 the three 
masses start at the positions marked by 1, 2,3 that are joined by the dashed line. They 
reach the positions marked by l',2',3' joined by the dotted triangle at time t = ff ■ 
The motion in the time interval [0, completely determines the rest of the motion 
that can be obtained by suitable reflections and permutations of particles. At the time 
t = ^ the three masses assume their positions at t = 0 up to a cyclic permutation 
(“choreography”). 


The coefficients '^n,m and Om can be calcnlated relatively rapidly dne to the linearity of 
the problem. Hence we have chosen a maximal trnncation order of n = 99. The results 
are again too complicated to be presented here and can be found in [5]. We only show 
the hrst few terms of the FT series for the lowest eigenvalue: 

“ = '-y-S+l5<®2'‘+') + 2l76 .(26) 

For yU = 0 the these coefficients coincide with those given in m 20.2.25 if = 2g is 
taken into account. Figure 0] shows the A, /i-dependence of a for a parameter region 
where we expect convergence of the FT series. As an example of a periodic solution we 
have chosen the values A = 2, /r = —0.3 and calculated the corresponding FT series as 
well as a numerical solution of 11231) with the same initial conditions and eigenvalue as the 
FT solution. The result is shown in the hgures [5] and [6] and demonstrates that the FT 
method satisfactorily works for this example. We conjecture that the periodic solution 
is unstable in this case and hence the distance to the FT approximation increases with 

X. 

4 . FT series for the figure eight solution 

The hgure eight solution, see hgure [71 is a stable plane periodic solution of the three 
body problem with equal masses that has been found numerically hrst in [3] and later 
been established by a rigorous existence proof |1]. It has attracted considerable interest 
i, 0, uni, mi, and has also been generalized to larger numbers of bodies [12] . One 
remarkable feature of it is that all three bodies follow the same curve but with a constant 
phase shift (“choreography”). This is similar to the property of a spin wave in a hnite 
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spatial structure, see |B]. Nevertheless, an analytical solution is not known. The Fourier 
components of this solution decay rapidly in order, hence a direct approximation by an 
FT series seems possible at hrst sight. Moreover, the choreographic character of the 
hgure eight solution can be easily implemented into the FT series. 


Let us hrst consider the equation of motion. We have three bodies with equal masses 
m, treated as point-like particles, that interact via gravitational forces. Anticipating 
that we consider periodic solutions with a period T we choose to = ^ as the unit 
of time and Xo = as the unit of length, where 7 denotes the gravitational 

constant. Passing to dimensionless quantities, denoting the position of the three bodies 

I Xi/(t) \ 

by r^(t) = '' \ = xfit) Fiyfit)^ i/ = 1,2,3 and time derivatives by a dot, we 

yfit) 


obtain the following equations of motion 


ra - ri rg - ri 

Tl — -j-pr + ]-f 

Ta-riF Ts-ril 


(27) 


and analogously for fa and fs where the corresponding equations are obtained by 
cyclic permutations of 1,2,3 (henceforwards denoted by cP). This equation of motion 
is invariant under the inhomogeneous Galilei group in two dimensions (6 parameters) 
and Kepler transformations (1 parameter). By a “Kepler transformation” we mean the 
multiplication of all lengths by and all times by /i^, where /i > 0 is some parameter. 
This transformation underlies the third Kepler law and is also utilized in the above 
introduction of dimensionless quantities. Utilizing Galilei invariance one chooses an 
inertial system such that the center of mass is always hxed to the origin: 


ri -f r2 -f rg = 0 , 


(28) 


and consequently 


f 1 -f f 1 -f f 1 = 0 . 


(29) 


Further constants of motion are the total energy and angular momentum that has the 
value L = 0 for the hgure eight solution. 

The numerical hgure eight solution suggests an FT series that starts with the terms 


xfit) = Al l A sint -|- ... (30) 

yfif) = i? 2,2 sin 2f -f ... . (31) 

However, this ansatz leads to the problem that the r. h. s. of (^Th is not diherentiable 
at A = 0 and hence the FT series does not exist. Attempts to weaken the strict cor¬ 
respondence between Fourier orders and powers of A avoid this problem but cause new 
ones. Another problem is that the FT series usually yields a one-parameter family of 
periodic solutions, as in the example of section | 2 l but such a non-trivial family of hgure 
eight solutions is not known. By “trivial” families of hgure eight solutions we mean the 
transformations induced by the invariance group of mentioned above. Hence the 
FT series method is apparently not suited to treat the problem at hand. 
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Yet we start a new attempt by considering relative coordinates and defining 

ri 2 = ri - r 2 = ri 2 exp(i (^ 12 ) and cP . (32) 

After some steps the new equation of motion is written as 

ri2 = ri2 012 - 2 rf2 + cos(<y9i2 - (fsi) + cos^ipu - (P 23 ), (33) 

012 = (-2 ri2 012 - sin((^i2 - (^ 31 ) - sin((^i2 - (^ 23 )) , (34) 


and the other components obtained by cP. This equation of motion is equivalent to fl27)l 
under the condition 

S = ri2 + r23 + r3i = 0 . (35) 

It can be shown that fl33|) and fl34D imply S = 0, hence S must be a constant of motion 
since it is a periodic function of time. But it could assume any value. If the particular 
value S = 0 is assumed for some solution of fl3^ and fl341) . then the definitions 


ri = - (ri 2 — r 3 i) and cP for r 2 and r 3 (36) 

o 

yields a solution of (|27j). 

Next we make the FT ansatz 

00 00 

ruit) = A2n,2mCos(4 7rnt), (37) 

n=0 m=n 

r23{t) = ri2{t - 27r/3), r3i(t) = ri2(t + 27r/3), (38) 

00 00 

(pi 2 {t) = A cos(t) + YY ^2m+l 2m+l COs(2 TT {2n + 1) t) (39) 

n=l m=n 

T23{t) = ri2{t - 27r/3), ip3i{t) = ri2{t + 27r/3) . (40) 


Thus ri 2 is assumed to be an even cos series and (fu to be an odd one. 9 ?ij(t) could be 
written more generally to contain a non-zero mean value. This would be equivalent to 
a constant rotation of the figure eight curve. We decided to set this mean value to zero 
and consequently obtain a figure eight curve standing upright, see below. Moreover, A 
is “concretized” by identifying it with the amplitude of the cost mode of (^ 12 , i. e. , 

1 

A = — / ipi 2 {t) cost dt. (41) 

Jo 

The numerical solution of [1] corresponds to a value of 

A = Ao ^ 2.07568 . (42) 

It will be in order to justify the above FT ansatz by referring only to general 
properties of the figure eight solution, not to its precise form according to the 
numerical solution obtainable from |1]. This would support the claim that the FT 
series approximation is, in principle, independent from the well-known numerical 
approximation. We first observe that rfi—t) = —rfit), hence xfit) and 1/3 (t) can be 
expanded into sin series. Anticipating that our ansatz leads to a figure eight curve 
standing upright, xfit) must be an even sin series and y{t) an odd one. This follows 
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from the symmetry properties at t j = |, namely xfin/2 — t) = —xfi7i/2 + t) and 
yfin/2 — f) = yfin/2 + f). By virtue of choreography, all yfif) are odd Fourier series 
and all xfif) even ones for i = 1, 2, 3. 

Since r3(0) = 0 = ri(0) + r2(0), c. f. fl28|) . it follows by the symmetry of the hgure 
eight solution and fl2^ that ri(0) = r2(0) = — |r3(0). Hence ri 2 (t) = ri(t) — r 2 (t) 
has a stationary value at t = 0, actually a local maximum. This makes it (albeit not 
necessary but) highly plausible to choose ri 2 {t) and </5i2(^) as cos series. Due to its 
dehnition ri 2 (t) will be an even cos series thereby justifying the ansatz fl57)) . Recall that 
all i = 1, 2, 3, and hence also yi 2 {t) are odd Fourier series. Since yi 2 = 'ri 2 sin (^12 
the obvious way to guarantee this is to choose y}i 2 {t) as an odd cos series. This completes 
the justihcaltion of fl37|l and fl3^ . 

Note further that the choreographic ansatz in (l38ll and (l4n|l can easily be eval¬ 
uated by means of the addition theorem for the cos function. This ansatz implies 
rij = Hqo + 0{X), ipij = 0{X) and hence the r. h. s. of (l33|) and (jSS]) will be differ¬ 
entiable at A = 0. It turns out that this FT ansatz yields unique solutions for the 
coefficients Anm and Bnm that can be calculated by computer algebraic means up to 
any desired order that is only limited by the available computer resources. 

But we have got another problem: The condition fl3^ is in general not compatible with 
flTr|) to (HOj) . This can be easily seen by considering the limit A —?• 0, where all three 
relative distances approach the value Hqo which means that the bodies form an equilat¬ 
eral triangle. Hence one of the relative angles (pij should approach a non-zero value, for 
example ^ and not 0. Another consequence of the violation of fl5SD is that the angular 
momentum is no longer conserved. 

Summarizing, we have obtained a 1-parameter family of periodic solutions of fl3^ and 
fl3T|l . but unfortunately of un-physical ones. Nevertheless, there is still a chance of 
analytically approximating the hgure eight solution. If the FT series is inserted, the 
quantity S becomes a function of A. It can be shown that Sj^(A) always vanishes since 
it is an odd Fourier series and a constant of motion. On the other hand, S 2 ,(A) will 
only vanish for a discrete set of zeroes A*. For these values, the FT series should give a 
physical solution of (ET]), especially the value A = Aq should correspond to the numerical 
solution of jl]. Hence we are left with the problem of determining the zeroes of Sa;(A) 
independently of the value A = Aq obtained from the literature. We will hrst illustrate 
the results of a low order truncation of the FT series. 


The coefficients of the FT series up the 6th order read 


AoO — 3^^^, y4o2 — • 
H22 — H24 — 


Aaa — 


8 

31/3 


, A46 — 


3^/3 ^ 3 X 31/3 

“ 256 ’ 

3^/3 ^ _ 5 X 31/3 

’ 1^’ “ 4096 ’ 

29 X 31/3 ^ 31/3 

, Am = 


49152 


8192’ 


(43) 

(44) 

(45) 
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-B33 — 



1 


B55 — 


9 

20480 ■ 


This yields the following analytical approximations via fl36|l : 


xsit) ^ 

yfit) ~ 


(A^ - 32A2 + 512) A2 sin(2t) (A^ + 192) A^ sin(4t) 
4096 31/6 49152 3 V 6 

A (3A^ — 32A^ + 256) sin(t) 3A^sin(5t) 

256 S'/® ^ 4096 3'/* ' 


(46) 


(47) 

(48) 


This result seems reasonable since it correctly renders xfif) as an even sin series and 
yfit) as an odd one. The coefficient of sin 3f must vanish since otherwise it would lead to 
a non-zero value of yfit) +y 2 {fi) +yfit) via choreography. In order to check the degree of 
approximation we have plotted the orbit corresponding to flTTl) and fl48|) together with 
the orbit obtained by an FT series of order 30 anticipating the result explained below. 
For the first curve the FT parameter A is chosen as the smallest real zero Ag ~ 2.09531 of 
S(A) evaluated up to 6th order in A. The coincidence between both curves is not perfect, 
see figure [HI but the example demonstrates that already an FT series approximation of 
low order gives a qualitatively correct picture of the figure eight solution. By the way, 
the main source of the deviation is not the restriction to the order 6 of the Fourier series 
but the error of the Fourier coefficients due to the relatively large value of Ag. 


Next we extend the calculations to the 30th order of the FT series. The detailed 
results for the series coefficients are too complicated to be presented here and can be 
found in [5]. But we will show the graph of the function Sa:(A), obtained by truncations 
of the FT series of order n = 8 to n = 30, see figure [9l This graphs makes it plausible 
that Sa;(A) has a real zero close to Ag, see fi42|) . corresponding to the physical value of 
the series parameter, but does not give a hint to the existence of further real zeroes. 
The zeroes A„ corresponding to truncations of order n can be numerically calculated 
and are shown in figure HOl It seems that they perform a damped oscillation about their 
limiting value for n —)■ cxo. Consequently, we take as an estimate of the true zero of 
Sa:(A) not Aq but the mean value of all A„, n = 12,14,..., 30 that amounts to 

Ai = 2.07552 , (49) 

and has a relative deviation from the value Ag obtained from [1] of less than 8 x 10“^. 
Our method to determine Ai could be criticized on grounds of its arbitrariness and 
crudeness but we think that it would be pointless to invoke more sophisticated methods 
that could only slightly improve the deviation from the “correct” value of the zero. 

The next task is to calculate the figure eight solution corresponding to the calculated 
series coefficients of order 30 and the estimate Ai of the physical value of the FT series 
parameter A. The result for the orbit has already be presented in figure [HI The difference 
to the figure eight curve obtained from the data of |1] would not be visible. Hence we 
have plotted the Euclidean distance a function of t, see 

figure [TTl where rf^(f) denotes the analytical approximation based on the FT series of 
order 30 and rf^(t) the numerical figure eight solution according to the initial values 
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Figure 8. Comparison of the figure eight solution curve corresponding to the FT 
series dni),® of order 6 (red dashed curve) with that of order 30 (blue curve). The 
FT parameter A of the first curve is chosen as A = Ae ~ 2.09531. Note that in contrast 
to figure 0 the figure eight is standing upright due to our choice of the constant value 
for the FT series of (pi 2 in (|5n)) . 
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Figure 9. Approximate graph of the function S 2 :(A) defined in (1551) corresponding to 
various orders n = 8,..., 30 of truncation of the FT series. The vanishing of Sa;(A) is a 
necessary condition for the FT series becoming a physically meaningful approximation 
of the figure eight solution. Accordingly we expect that the FT series parameter A 
should assume a value close to 2.07552, see also figure [TUI 
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Figure 10. Numerical zeroes A„ of Sa;(A) based on truncations of the FT series of 
order n = 12,..., 30. The mean value Ai = 2.07552 gives an estimate of the correct 
zero. 


published in [1] . It is clear that we had to modify these initial values by a scale (“Kepler”) 
transformation and by suitable reflections for the sake of comparison with the FT series 
approximation. The result is a maximal deviation of 1.6 x 10“^ between both solutions 
that is larger than expected but explains why the difference between both hgure eight 
curves would not be visible. 
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Figure 11. Plot of the Euclidean distance as a function of time t between the analytical 
approximation based on an FT series of order 30 and the numerical figure eight solution 
based on the initial values published in [3]. 


5. Summary and Outlook 

We have proposed a method suited to analytically approximate periodic solutions of 
non-linear equations of motion, the Fourier Taylor (FT) series. If the method works it 
yields the series coefficients in a recursive way such that each step involves a new Fourier 
order and simultaneously improves the accuracy of the old Fourier coefficients. Typical 
applications of this method include mechanical oscillations about some ground state and 
extend the validity of the approximation beyond the linear regime, as in the example of 
section [2j Even for linear problems the FT method can be used as a simpler alternative 
to perturbation theory, see section [3l In these cases the FT series yields a family of 
solutions depending on a parameter A and the domain of application is limited by the 
convergence radius w. r. t. A and, additionally, by the possible truncation order of the FT 
series. In most cases the FT coefficients can only be calculated by the aid of computer 
algebra and hence the maximal truncation order is given by the available computer 
sources. Nevertheless, the FT series approximation has some advantages compared 
with a purely numerical treatment of the problem: For example, the FT coefficients 
may depend on further parameters of the problem and thus an FT approximation may 
be equivalent to an inhnite number of numerical calculations for different values of the 
parameter, just as in the case where a closed formula for a class of solutions has been 
found. 

From this point of view the application of the FT series method to the hgure eight 
problem is somewhat untypical. We could not directly apply the method to the equation 
of motion but only to another equation for the relative coordinates that follows from 
the original one. Conversely, the new equation implies the original one only if some 
addition condition (l3^ is satished. The FT series violates this condition, and thus 














FT series for the figure eight solution 


17 


produces un-physical solutions, except for a particular value Ai of the series parameter 
A. For this value the FT series approximates the hgure eight solution known from |1], 
where the error depends on the truncation order. In principle it can be arbitrarily small, 
in practice we have to be content with an error of 1.6 x 10“^ for a truncation order of 
n = 30. 

One may ask: What is the particular virtue of the FT series in the hgure eight example, 
except for illustrating the scope of the method? There exist already two ways to access 
the problem, namely the numerical solution found by some educated guess Q. Q and 
the mathematical existence proof |1]. Compared with both ways the FT series method 
seems to be inferior: The result obtainable in practice is less accurate than the numerical 
one and it is not rigorous in so far as the convergence of the FT series has not been 
proven. We think, however, that the virtue of the FT series method is, hrst, that it 
provides a third approach to the hgure eight problem that is, in principle, independent of 
the two other ones. For this reason we have carefully identihed the assumptions leading 
to the FT ansatz in section 01 Secondly, the FT series method requires less ingenuity, 
so to speak, compared with the two other methods. Once the adequate ansatz has been 
found, the method works automatically and can (and should) be performed by computer 
software. Hence in this way it might be possible to easier hnd new solutions of the three 
body problem that are conjectured to exist and to have special symmetries. But this is 
a task for future work and beyond the scope of the present paper. 
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